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ABSTRACT 

We investigated the absolute age of the Galactic globular cluster M71 (NGC 6838) by using optical ground- 
based images {u',g',r',i',z') collected with the MegaCam camera at the Canada-Erance-Hawaii-Telescope 
(CEHT). We performed a robust selection of field and cluster stars by applying a new method based on the 3D 
(r', u' - g', g' - r') Color-Color-Magnitude-Diagram. The comparison between the Color-Magnitude-Diagram 
of the candidate cluster stars and a new set of isochrones, at the locus of the Main Sequence Turn Off (MSTO), 
suggests an absolute age of 12+2 Gyr. The absolute age was also estimated using the difference in magnitude 
between the MSTO and the so-called main sequence knee, a well defined bending occurring in the lower main 
sequence. This feature was originally detected in the near-infrared (NIR) bands and explained as a consequence 
of an opacity m echanism (collisio nally induced absorption of molecular hydrogen) in the atmosphere of cool 
low-mass stars (iBono et al. Il201(]l) . The same feature was also detected in the r' , u' - g' and in the r',g' - r' 
CMD, thus supporting previous theoretical predictions bv iBorvsow et al. I (119971) . The key advantage in using 
the as an age diagnostic is that it is independent of uncertainties affecting the distance, the reddening and 
the photometric zero-point. We found an absolute age of 12+1 Gyr that agrees, within the errors, with similar 
age estimates, but the uncertainty is on average a factor of two smaller. We also found that the is more 
sensitive to the metallicity than the MSTO, but the dependence becomes vanishing using the difference in color 
between the MSK and the MSTO. 

Subject headings: globular clusters: individual (M71) 


1. INTRODUCTION 

The Galactic Glo bular Cluster (GGCs) M71 -NGC 6838- 
([Fe/H]=-0.78 dex. lHarrIs1ll996h is a very interesting stel¬ 
lar system since it belongs to th e small sampl e of metal- 
rich GCs pres ent in our Galaxy ([Harris 1 1199^ and in the 
Local Group (ICezario et al. 1 12013lk Moreover, it is lo¬ 
cated relative ly close, and indeed its distance is smaller 
than 4 Kpc (Grundahl, Stetson & Andersen 1 120021) . This 
means that M7 lisa fundamental laboratory t o constrain the 
evolu tion of old, metal-rich, low-mass stars (iHodder et al. I 
119921) . The absolute age of M71 can also play a key role 
in the occurrence of a bifurcation in the age-metallicity re- 
lation of inner and outer halo GGCs recently suggested by 
Potter, Saraiedini & Anderson] (1201 Ih . Moreoyer, M71 to- 
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gether with a few more metal-rich GGCs (47 Tuc, NGC 6528, 
NGC 6553) plays a fundamental role in the calibration of 
spectrophotometric indices adopted to constrain the age, the 
metallicity and the shape of the Initial-Mass-Eu nction of unre¬ 
solved stellar populations in e a rly type galaxies dBoselli elalJ 
20091 ICappellari et al. II2012I IConrov & van Dokkum (20121: 


Conroy, van Dokkum & Graves l2013USpiniello et al. rioi2h . 

The main drawback of M71 is that it is affected 
by a large reddening wi th a mean valuJ^ of < 
EiB- V) >=0.323+0.017 mag (ISchlegel, Einkbeiner & Davis I 
119981) and by a differential reddening with a mean value 
of < d E(B-V) >=0.035+0.015 mag, which extents 
up to the maximum value of 6 E(B-V)-Q.Q74 mag 
( iBonatto, Campos & Oliveii^ 120131) (for the red- 


1976; 


dening values see al so iKron &Yjuetter| 

Frogel, Persson, & Cohen I 119791 [ktodder et al. I _ 

Kraft & Ivans Il2003h . This is mainly due to the fact that M71 


19921 


is located at very low Galactic latitude (1=56°.75, b=-4°.56). 
The literature concerning the age estimates of M71 is quite 


rich. By comparing optic al CCD photometry with the oxy ¬ 
gen enhanced mode ls of Bergbush & VandenBergI (119921) . 
[Hodder et al. I (Il992h found that M71 has an age in good 
agreement with the 14+2 and 16+ 2 Gyr isochrones. They 
also found that M71 is c oeval with 47 Tuc, while b y using the 
same theoretical models iGeffert & Maintz I (l2000h suggested 
for M71 an older age (18 Gyr). By calibrating the cluster 
distance with the Hipparcos field subdwarfs p arallax e s, and 
by adopting a reddening of E(B-V)-0.2H mag. iReidI (119981) 
found a true distance modulus of ju=13.19+0.15 mag and an 
absolute age of 8+1 Gyr (as estimated from his Fig. 10a). He 
also suggested that 47 Tuc, within the uncertainties, is coeval 


The mean value has been calculated by the online tool available at 
http://irsa.ipac.caltech.edu/applications/DUST/ 
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(10+1 Gyr). 

The Hipp arcos catalog was also used with Stroemgren pho¬ 
tometry by Grundahl. Stetson & Andersen I (l2002h . who de¬ 
rived a true distance modulus of ju=12.84+0.04+0. iQ mag 
by adopting E(B-V)-0.2 % mag. By comparing the M 71 data 
with the isochrones from IVandenBerg et al.1 (120001) . the au¬ 
thors also provided an absolute age of 12 Gyr and found the 
same age for 47 Tuc. 

After the transformation of the IClem et ^ (l2008h M71 
fiducial sequence t o the SPSS bands with the equa¬ 
tions provided bv jTuokeret alJ (l2006l) . and assuming 
E(B-V)-0.2S mag, lAn et al. I (120091) found a distance mod¬ 
ulus of yu= 12.86+0.08 mag and 12.96+0.08 mag re- 
spectively according to the metallicity scale prov ided by 
iKraft & Ivans I (120031) and ICarretta & Gratton I (119971) . 

To fu rther improve the temp erature sensitivity around the 
MSTO, iBrasseur et al. I (1201 (il) adopted optical-NIR colors 
and, by using a true distance modulus of //=13.16 mag and 
a mean cluster reddening of E(B-V)=Q.2Q mag, they found 
an absolute age of 11 Gyr. This age estim ate was based on a 


new set of cluster isochrones computed by 

VandenBerff et al. 1 

(I2012^. More recently. IVandenBers et al. 

(|2013l hereinafter 


VI3) by using homogeneous optical photometry based on 
WFC/ACS images, provided accurate relative ages for 55 
Galactic globulars. They adopted an improved vertical 
method and a new approach to fit Horizontal Branch (HB) 
stars with predicted Zero Age Horizontal Branch models. 
They found for M71 an absolute age of 11.00+0.38 Gyr. 

In a more recent paper, by inve s tigatin g the age of 61 GCs, 
Teaman, VandenBerg & Mendel I (1201 3L hereinafter LVM13) 
found that for abundances [Fe/H]>1.8 dex, the GCs obey to 
two different Age-Metallicity-Relations (AMRs). In particu¬ 
lar, they found that one-third of the entire sample is, at fixed 
cluster age, systematically more metal-rich by 0.6 dex. More¬ 
over, they suggest that the bulk of the metal-rich sequence 
formed in situ in the Galactic disk, while a significant frac¬ 
tion of the metal-poor globulars formed in dwarf galaxies 
that were subsequently accreted by the Milky Way. However, 
the current scenario concerning the absolute age distribution 
of GCs is far from being settled. In a similar investigation 
iMarin-Franch et al~l (l2009t) . by using homogeneous photom¬ 
etry for a sample of 64 globulars, found that the bulk of GCs 
are coeval, the spread in age being smaller than 5%, and do 
not display an AMR. They also found a small group of sys¬ 
tematically younger clusters that display an AMR similar to 
the globulars associated to the Sagittarius dwarf galaxy. To 
take account of the observed age distribution, they suggested 
that the old globulars formed in situ, while the young ones 
formed in dwarf galaxies that were later accreted. 

V13 and LVM13 did not include in their sample red¬ 
dened bulge GCs, but on the basis of the AMR relation they 
found for metal-rich clusters they predicted steadily decreas¬ 
ing ages for more metal-rich bulge clusters (see their sec¬ 
tion 7.3). In particular, for the prototype metal-rich bulge 
cluster NG C 6528 they predicted an age of 10 Gyr. On the 
other hand, iLagioia et al. I (l2014t) by using WFC/ACS and 
UVIS/WFC3 images that cover a time interval of 10 years 
performed an accurate proper motion separation between field 
and cluster stars in NGC 6528. On the basis of optical pho¬ 
tometry they found that NGC 6528 is old (12 Gyr) and coeval 
with 47 Tuc. 

We used t(ie ICardelli. Clavton & Mattiis 1 119891) law to unredden the 
original value of My=\3Jl mag. 


The above evidence indicates that we still lack solid con¬ 
straints on the age distribution of metal-rich GCs. In this in¬ 
vestigation we address the absolute age of M71 by adopting 
optical ground-based (SDSS bands) collected with MegaCam 
at Canadian-France-Hawaii-Telescope (CFHT). The layout 
of the paper is the following. 0 In section we present the 
dataset and the strategy adopted for data reduction. In section 
[3 we discuss the approach we adopted to identify candidate 
cluster and field stars. The theoretical framework adopted to 
estimate the absolute age of M71 is presented in section|4]and 
12 while in section |2 we discuss the absolute age estimates 
based on the classical main sequence turn-off and on the main 
sequence knee. Finally, in section |7] we summarize the results 
and outline future perspectives. 

2. OBSERVATIONS AND DATA REDUCTION 

To investigate the absolute age of M71 we used the images 
of the MegaCam (36 Charge Coupled Devices, CCD; total 
field of view (FoV);l°xl°; scale: 0.187" /pixel) mounted on 
the CFHT. The dataset includes 50 dithered imageJ^ taken 
with the u'g'r'i'z' (SDSS) bands (iFukugita et all 19^ . We 
retrieved from the Canadian Astronomy Data Centre, five 
shallow and five deep images per filter with the following 
exposure times (ETs): ET(u')—30, 500 sec; ET(g',r')—5, 
250 sec; ET{i')^5, 300 sec; and ET(z')^15, 500 sec. The 
mean seeing ranges from ~1.0 (u',g',z') up to ~1.3 ()',/) 
arcsecond s. The data were pre-process ed with the ELIXIR 
programs (iMagnier & Cuillandre 1 120041) . The globular M71 
covers a small sky area, and indeed its tidal radius is r,-%'.9 
(iHarris I Il996h . therefore, we only used the 20 innermost 
CCDs for a total EoV of 0°.6xl°.0. Eor each frame we pro¬ 
vided an accurate PSE (Point Spread Euncti on) photometr y 
by using the DAOPHOT IV and ALLSTAR (iStetson 11198% . 
and we used DAOMATCH/DAOMASTER to scale individual 
chips on a common geometrical reference system. Once we 
obtained the global catalog ( master list of sta rs) of the entire 
dataset we run ALLERAME (lStetsonlll994h simultaneously 
over the entire set of images. We obtained a list of 370,000 
stars with at least one measurement in four different photo¬ 
metric bands. To calibrate the instrumental mag nitudes, we 
used th e local standard stars (~6,000) provided bv IClem et alJ 
(|2003). To validate the adopted transformations, EigUshows 
the comparison in four different Color-Magnitude-Diagrams 
(CMDs) between the curren t photometry and the ridgelines 
(red solid lines) provided bv IClem et al.l (Hool). They agree 
quite well not only along the Red Giant Branch (RGB), but 
also along the Main Sequence (MS). 

The anonymous referee noted that for r'-band magnitudes 
fainter than ~22 mag {r',g' - r' CMD; r' ~20.5 mag for 
the r', u' - g' CMD) the ridgelines provided by Clem attain 
colors that are slightly bluer (redder for the r', u' - g' CMD) 
than observed. The difference is mainly caused by the differ¬ 
ent approaches adopted to perform the photometry, in select¬ 
ing candidate cluster stars and in the calibration to standards. 
Our approach is based on simultaneous photometry on shal¬ 
low and deep u',g', r', i',z' images. This allowed us to reach 
a better photometric precision in the faint magnitude limit. 
Stars plotted in Eig. 1 have been selected according to the 
cluster radial distance (60<r<250 arcsec), and to photometric 
quality parameters (cr^'-r- = 0.07 mag, x ~ 1, \sharp\ < 1). 

The cun'ent optical images were acquired during two nights: July 8 and 
13, 2004; proposal IDs: 04AC03, 03AC16; P.I. J. Clem; observing program: 
’CFHT Star Cluster Survey’. 
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Data plotted in this hgure indicate that the lower MS is bet¬ 
ter defined, less affected by the field star contamination (see 
Fig. 12 in Clem et al. 2 008). We adopted the same local 
standard stars prov ided bv IClem et al.l (l2007h and followed a 
similar procedure (IClem et al.ll2008h to transform the instru¬ 
mental magnitudes into standard magnitudes. Unfortunately, 
we cannot perform a detailed comparison among the different 
calibration equations, since their zero-points and coefficients 
of the color terms are not available. 

The intrinsic error in magnitude and in color are plotted as 
red error bars. They take account of the photometric error 
and of the absolute calibration and attain values of the order 
of a few hundredths of magnitude down to the lower main 
sequence. This trend is also supported by the photometric 
error in the r'-band plotted in the rightmost panel of the same 
figure. It is of the order of a few hundredths of magnitude 
down to the bending of the MS (r' ~23 mag). 

Finally, we performed an astrometric solution of the catalog 
by using the UCAC4 catalog. The root mean square error of 
the positions is cr ~55 10“^ arcseconds. 

3. IDENTIFICATION OF FIFED AND CLUSTER STARS 

To separate candidate cluster and field stars we devised a 
new approach that fully exploit the multi-band data set we are 
dealing with. We performed a preliminary radial selection and 
to avoid the crowded central regions we only selected stars lo¬ 
cated inside an annulus with r=2'and r=5'centred on the clus¬ 
ter. On the basis of this sample we generated two isodensity 
maps in the /, u' - g' and in the r',g' - r' CMDs (see Fig.|2]l. 
To properly constrain the cluster ridge line we developed a nu¬ 
merical algorithm that pin points the peaks of the isocontour 
plots and provides a preliminary version of the cluster ridge¬ 
line. Note that in the approach we devised, the overdensity 
caused by red HB stars was neglected. The preliminary ridge¬ 
line is then fit with a bicubic spline and visually smoothed, 
in particular in the bright portion of the RGB. The analytical 
ridgelines are sampled at the same r' magnitude levels (see 
TablellJ. The red lines plotted in the left and in the right panel 
of Fig.|2]show the final version of the estimated ridgelines. 

The two analytical ridgelines were used to generate a 
3D plot Color-Color-Magnitude Diagram (CCMD)- r',u' - 
g',g' - r' -of the cluster. The Fig. [3] shows the 3D ridgeline 
(red solid line) together with the entire sample of stars with 
measurement in the the u', g', r' bands. We considered as can¬ 
didate cluster stars those located within the 3D ridgeline ±\cr, 
where cr is defined as two times the quadratic sum of the pho¬ 
tometric errors in the three quoted bands. 

The above approach has several indisputable advantages 
when compared with classical photometric methods in select¬ 
ing field and cluster stars. 

i) The selection in the CCMD takes advantage of the typi¬ 
cal effective temperature correlation of stellar structures in a 
color-color plane. Current approach takes also advantage of 
the fact that we are simultaneously using photometric bands 
covering a broad range in central wavelenghts, and in partic¬ 
ular, of the w'-band. This means a solid separation not only 
between sta rs and field galaxies , but also between field and 
cluster stars iBono et al. I (I2010bl) . However, the color-color- 
diagram approach is prone to degeneracy between dwarf held 
and giant cluster stars sharing very similar optical colors. 

ii) The selection in the CMD takes advantage of the typi¬ 
cal correlation between cluster stars located at the same dis¬ 
tance (apparent magnitude) and their color. However, the 
CMD approach is prone to spurious selections interlopers lo¬ 


cated at different distances, i.e. held stars associated either 
to the Galactic halo or to the Galactic disk stars. This means 
that after the selection of the stars located across the ridge¬ 
line we need to take account of the fraction of held stars that 
have been erroneously classihed as candidate cluster stars (see 
iPi Cecco et al. lIMo 

The approach based on the CCMD takes account of the ad¬ 
vantages of both the color-color-diagram and of the CMD se¬ 
lection criteria. To validate the approach we adopted. Fig. |4] 
shows the r',g' — r' CMD, before (left) and after (middle) the 
selection. The candidate held stars are plotted in the right 
panel. The magnitude and color distribution of held stars sup¬ 
port the plausibility of the adopted selection criterion, since 
we recover the two typical pea ks in g' - r' color at ~0.6 and at 
~1.5 mag of held dwarf stars (llvezic et al. 1120081) . Data plot¬ 
ted in the middle panel display a smooth and well sampled 
CMD from the tip of the RGB down to the lower MS. More¬ 
over, the presence of held stars in the same region of the CMD 
in which we would expect candidate cluster stars on the basis 
of the CMD, is further supporting the approach we developed. 

We highlight that the isodensity maps plotted in Fig. 2 
show a secondary peak at magnitude slightly brighter (17< 
r' <17.5 mag) and bluer («' - g' ~\.2 mag; g' -r' ~ 0.6 mag) 
than the ridgeline of the candidate cluster stars. We performed 
several tests changing the ridgeline either in the r' ,u' - g' or 
in the r' ,g' - r' density map and checking the impact on the 
CCMD. We found that this blue plume mainly includes can¬ 
didate field stars. This finding is supported by the smooth dis¬ 
tribution of cluster MSTO stars in Fig. 1 and by the presence 
of the blue plume stars in the CMD of candidate field stars 
plotted in the left panel of Fig. 4. We thank the anonymous 
referee for drawing our attention on this group of stars. 

4. THEORETICAL FRAMEWORK AND CLUSTER ISOCHRONES 

To estimate the cluster age we used the evolutionary 
tracks and the cluster isochrones of the Pisa Stellar Evo¬ 
lution Data Bas43 for low-mass st ars, computed using 
FRANEC stellar evolutio nary code (iDeglTnnocenti et al. I 
120081 iTbgnelli et ffi~ll201 Ih . The input physics and the phys¬ 
ical assumptions adopted to c onstruct the evolutionary track s 
have already been discussed in iDeirOmodarme et al. I (l2012h . 
while an analysis of the main theoretical uncertai nties affect¬ 
ing stel lar tracks and isochrones is provided in IValle et alTI 
(I2013al[bl) . Among the available grids of models, we take 
account of those computed with an a-enhanced chemical 
mixture ([g /Eel--F0.3) and the re cent solar heavy-element 
mixture by lAsplund et al. I (120091 As09). Recent spectro¬ 
scopic estimate s give for M71 an iron abundance of [Ee/H]=- 
0.82+0.02 dex (ICarretta et al. Il200^ Ca09). This abundance 
is based on a reference solar iron abundance of log eplEel) = 
7.54. To transform it into the lAsplund et al. I (120091) solar iron 
abundance (log^(Ee)=7.50), we used the following relation¬ 
ship; 

[Ee/H]As 09 = [Fe/H]ca 09 + log eo(Ee)ca 09 - log eo(Fe)As 09 = 

= [Fe/H]ca 09 + 0.04 (1) 

This means an iron abundance for M71 of [Fe/H]=-0.78 
dex, and by assuming an a-enhancement of [a/Ee] = +0.3 
dex, a global metal abundance per unit mass of Z=0.0037 
and a primordial helium abundance of 7=0.256. The clos¬ 
est metallicities for a-enhanced isochrones currently avail¬ 
able in the Data Base are: Z=0.002, 7=0.252; Z=0.003, 

http://astro.df.unipi.it/stellar-models/ 
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y=0.254 and Z=0.004, 7=0.256. The helium abundances 
as a function of the metal abundances were estimated by us- 
ing a linear helium-to-metal enrichment ratio of AT/AZ = 2 
(iFukugita & Kawasaki 11200^ . The above chemical composi¬ 
tions imply iron abundances of [Fe/H]=-1.04, [Fe/FI]=-0.87, 
and [Fe/H]=-0.74 dex, respectively. The adopted sets of evo¬ 
lutionary models bracket, within the errors, the observed iron 
abundance of M71. 

The two adopted grids of evolutionary models cover the 
typical range in mass of low-mass stars, namely M- 0.SO¬ 
LI 0 Mq. Moreover, they were constructed by adopting three 
different values of the mixing length parameter; ml-l.l, 1.8, 
and 1.9. The corresponding cluster isochrones cover the age 
range from 8 to 15 Gyr. The luminosities and the effec¬ 
tive temperatures provided by the evolutionary models were 
transformed into ma gnitudes and color indexes by using the 
synthetic spectra bv IBrott & Hauschildt I (I2005h for 2000 K 
< Tg/y < 10000 K and by iCastelli & Kuru(Sl (I2003h for 
lOOOOK < Teff < 50000 K. The transformations into 
the observational pl ane follow the prescriptions discussed in 
iGirardi et al I (1200^ . 

To compute the AB magnitudes in the u'g'r'i'z' photomet¬ 
ric system we used the USNO40 Response Function^ with 
an air mass of 1.6. This is the air mass value at which the 
observations of th e local standard stars in M71 provided by 
IClem et aH (l2007h were performed. 

Fig. 0 shows the comparison between candidate cluster 
stars and the set of isochrones constructed assuming by 
Z=0.003, 7=0.254 and a ml-1.9. The cluster isochrones were 
plotted by adopting a true distance modulus of /t= 13.07 mag 
and a cluster reddening of E{B-V)-0.25 mag. The adopted 
values agree quite well with similar estimates available in 
the literature (see Table |2]i. The selective absorptions in 
the individual bands were estimated by assuming a ratio 
of absolute to selective extinction of Rv=3.1 and the em - 
pirical reddening law bv ICardelli, Clayton & Mathis I (Il989h . 
The extinction ratios were computed by adopting the effec¬ 
tive wave kngths of the u'g 'r'i'z' photometric system pro¬ 
vided by (ISmith et al. Il2002h . For the five adopted bands we 
found the following selective absorption ratios: A„//Av=1.58, 
A„-/Av=l .20, A.>/Av=0.87, A,v /Ai/=0.66. and A-/Av=0.49 
(see also iDi Cecco et aT~l 1201 Oh . The black arrows plotted 
in the bottom left corners display the reddening vectors in 
the four different CMDs. The two red lines display the clus¬ 
ter isochrones for 11 and 13 Gyr. Theory and observations 
agree quite well not only along the evolved sequences (RGB; 
Sub Giant Branch, SGB), but also along the MS. There is a 
mild evidence that cluster isochrones in the lower main se¬ 
quence (r' <20.5 mag) become slightly bluer than observed 
stars in the r', g' - r' (panel b), r', g' - i' (panel c) and 
in the r', g' - z! (panel d), but the agreement across the 
turn-off regions is quite good in the four different CMDs. 
Similar discrepancies between observed and predicted col- 
ors have also been foun d in the optical and in the NIR bands 
jKucinskas et al. l[2005b . The only exception is the r', u'-g' 
CMD (panel a), in which the observed stars in the MSTO 
region (r' ~17.8 mag) are a few hundreds of magnitude 
blu er than pr e dicted . A similar discrepancy was also found 
by lAn et al. I (|200^ using g-enhanced cluster isochrones 
computed with YREC code dSills, Pinsonneault & Terndrupj 

They are available at the following url 
http ://www-star.fnal.go v/ugriz/Filters/response.html see for more details 
(see lSmith et al. 120021) 


l2000t iDelahave & Pinsonneault I [20061) and transformed into 
the observatio nal plane by adopting M ARCS stellar atmo¬ 
sphere models (iGustafsson et al. Il2008h . They suggested that 
the main culprit could be missing opacities in the model at¬ 
mospheres at short wavelengths, as well as the color trans¬ 
formation for the highly reddened cluster M71. Current find¬ 
ings support the discrepancy, but the comparison also display 
a good agreement along the MS (r' <20 mag). Thus suggest¬ 
ing that the mismatch between theory and observations might 
also be due to marginal changes in the adopted chemical com¬ 
position. 

To further constrain the above working hypothesis, we per¬ 
formed the same comparison, but by adopting a slightly more 
metal-rich chemical composition (Z=0.004, 7=0.256). Data 
plotted in Fig. |6] clearly show that the discrepancy between 
theory and observations increases. This applies not only to the 
MSTO region where the difference is of the order of 0.3 mag, 
but also along the RGB/SGB and the MS, and indeed the 
isochrones attain colors that are systematically bluer than ob¬ 
served. Moreover, it is worth mentioning that the increase of 
0.13 dex in iron abundance also causes a clear discrepancy 
along the SGB region in the other CMDs (panels h, c, d) and 
does not improve the agreement in the lower MS. 

Current evidence further supports the strong sensitivity of 
the u'-band to the metal content and suggests that the mild 
discrepancy found in the color of TO region of the r', u' - g' 
CMD might also be affected by the adopted abundance of a- 
elements and/or of CNO-elements. 

5. CONSTRAINING THE NATURE OF THE KNEE IN THE LOWER 

MS 

The occurrence of a well defined knee in the low-mass 
regime of the MS has already been brought forward in the 
literature. It has already been detected in several old (cu Cen, 
Pulone et al. 1998; M4, Pulone et al. 1999, Milone et al. 
2014; NGC 3201, Bono et al. 2010; 47 Tuc, Lagioia et al. 
2014; NGC 2808, Milone et al. 2012) and intermediate-age 
dSaraiedini, Potter & Kirkpatrick I l^09h stellar system s and 
in the Galactic bulge (IZoccali, Cassisi & Frogeni2000h by us¬ 
ing NIR CMDs. Current evolutionary prescriptions indicate 
that this feature is mainly caused by th e collisionally induce d 
absorption of H 2 at NIR wavelengths (iSaumon et al. |[l994h . 
The shape of the bending marginally depends on the metal 
content, but the magnitude of the knee is essentially indepen¬ 
dent of cluster age and of metallicity. 

The above theoretical and empirical evidence suggests that 
the magnitude difference between the MS knee (MSK) and 
the MSTO, i.e. is a robust diagnostic to constrain the 

absolute cluster age. The key advantages of the new approach 
is that the above difference in magnitude is independent of 
uncertainties affecting both the distance and the reddening of 
the stellar system. Recent empirical evidence indicate that the 
this new approach provides cluster ages that are a factor of two 
more pr ecise when compared with the classical method of the 
MSTO (iBono et al. II 2 OIOI; Sarajedini, Potter & Kirkpatrick I 
1200 ^ . 

In passing, we also note that the method appears more 
robust than the traditional horizontal (difference in color be¬ 
tween the RGB, typically 2.5 mag brighter than the MSTO, 
and the MSTO ) and vertical (difference in magnitude be¬ 
tween the HB luminosity level, typical ly at the RR Lyrae in ¬ 
stability strip, and the MS TOl methods ([S tetson et al. f (ll999l) : 
[Buonanno et al. I (Il998h : iMarin-Franch et al. I (120091) : V13). 
According to theory, MS stellar structures with a stellar mass 
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of M ~0.5-0.4 Mq are minimally affected by uncertainties 
in the treatment of co nvection, since the convec tive motions 
are nearly adiabatic (ISaumon & Marled l2008l) . The same 
outcome does not apply to cool HB stars and to SGB stars 
adopted in the vertical and in the horizontal method, respec¬ 
tively. Moreover, the MSK can be easily identified in all stel¬ 
lar systems with a well populated MS. It is independent of 
the uncertainties affecting the estimate of the HB luminos¬ 
ity level when moving from metal-poor (blue and extreme 
HB morphology) to m e tal-rich (red HB morp hology) GCs 
(iCalamida et al. I (I2007h : Hannicola et al. I (I2009h : V13). Mo- 
rover, the method is also independent of the theoretical 
uncertainties plaguing the color-temperature transformations 
required in the horizontal method to constrain the cluster age. 

In this context it is worth mentioning that both the verti¬ 
cal and the horizontal method are robust diagnostics of op¬ 
tical CMDs. Recent accurate and deep NIR CMDs show 
that HB stars display a well defined slope when moving from 


hot and extreme HB stars to cool red HB stars (|Del Principe 

2006 

:ICoDDola et al. 1201 lUMilone et al. l2013HStetson et al. 

2014 

). The same outcome applies to the CMDs based on near 


UV and far UV bands (iFerraro et al. 1120 12h . This means that 
the identification of the HB luminosity level required by the 
vertical method is hampered by the exact location in color of 
the anchor along the HB. Moreover, the difference in color 
between the MSTO and an anchor along the RGB, required 
by the horizontal method, is hampered by the fact that the 
RGB becomes almost vertical in NIR bands. This means that 
the difference in color is steadily decreasing (iCoppo l a et al. I 
IMTllStetson et al. lIMIl) 

Finally, the minimal dependence on the metal content allow 
us to tight correlate the directly to the absolute age of 
the stellar system. This means that the new approach provides 
absolute age estimates, while the horizontal and the vertical 
method do provide estimates of the relative age. 

5.1. The impact of CIA on the knee of the lower MS 

To further constrain, on an empirical basis, the robustness 
of the to estimate the absolute age of GCs we also in- 
vestigated the beding in the optical bands (iBono et al. 1120 10^ 
and in particular in the u'- and in the g'-band. The reason 
is twofold. 1) Detailed atmosphere models taking account of 
the collision-induced absorption (CIA) opacities of both H 2 - 
H 2 and H 2 -He indicate that CIA induces a st rong continuum 
depression in the NIR region of the spectrum (iBorvsow et al. I 
119971 [Borysow, Borvsow & Fu II 2 OOOI: iBorvsow II2002I) . The 
same outcome applies to the UV region of the spectrum for 
surface gravities and effective te mperatures typical of cool 
MS stars (see Figs. 5 and 6 in IBorvsow et al. Ill997l) . It 
is interesting to note that the above effect is also anticorre¬ 
lated with the iron abundance. The increase in the iron abun¬ 
dance causes, at low effective te mperatures, an increase of the 
molecular opacities (see Fig. 7 in iBorvsow et al. I[l997h . lead¬ 
ing to a weaker impact of the CIA on the continuum. 2) Inter¬ 
estingly enough, recent deep and accurate optic al-UV CMD 
already showed the presence of a well define d knee in the 
low-mass regime of the MS. I An et al. 1 (l2008h clearly iden¬ 
tified a well defined knee in the lower MS of two old open 
cluster; NGC 2420 and M67 (see their Fig. 16). A similar evi¬ 
dence was also found in the V,B-I CMD of M4 (iStetson et al. I 
|2014|) . The lack of a clear evidence of a knee in the lower MS 
can be associated either to the use of V,R,I bands in which the 
knee is less evident or to the fact that the photometry is not 
deep/accurate enough to properly identify the knee. 


To properly constrain the impact that CIA has on the UV 
flux, and in particular on the blue magnitude and colors, we 
performed a specific test. We adopted an old (12 Gyr) clus¬ 
ter isochrone and selected three stellar structures with stellar 
masses ranging from 0.30 to 0.45 Mq. The surface gravities 
and the effective temperature of the selected models are listed 
in Table |3] Using these values as input parameters, we com¬ 
puted two specific sets of atmospheric models (spectra), one 
that neglects and one that takes account of the CIA opacities. 

Data plotted in Fig. Qdisplay that CIA has a substantial im¬ 
pact on the emerging flux of high-gravity, low effective tem¬ 
perature stellar structure. 

In the color-magnitude region common to the two sets, the 
models that take account of CIA (see red line in Fig. |7]i are, 
at fixed stellar mass, ~0.1 mag fainter than those without CIA 
opacity, the differences increasing as the effective temperature 
decreases. Notice also that, at a fixed luminosity, the inclusion 
of CIA produces models ~0.05 mag bluer than the models 
without CIA. The difference becomes more evident in the very 
low-mass regime and at longer wavelengths. 

In this context it is worth mentioning that the above at¬ 
mosphere models either taking account of or neglecting the 
CIA opacity rely on the input physics adopted in the lat- 
est version of the P HOENIX BT-Settl atmosphere structures 
(I Allard et al. 11201 lb . These models when compared with pre¬ 
vious atmosphere models computed by iBrott & Hauschildt I 
(120051) present two differences: a) the former use the solar 
mixture provided by Asplund et al. (2009), while the latter 
use the solar mixture provided by Grevesse & Noels (1993); 
b) the former models include a more complete list of CIA 
opacities. In particular, they take account of molecular hy¬ 
drogen (H 2 ), molecular nitrogen (N 2 ), methane and carbon 
dioxide (CO 2 ). The main difference with previous models is 
in the use of opacities provided by Abel & Frommhold (2011) 
that provide less absorptions at high temperatures when com¬ 
pared with earlier computations by iFu, Borvsow & Moraldi! 
(Il996h . A more detailed list of the adopted CIA opacities used 
in the current BT-Settl atmosphere models is given in Table |4] 
together with their references. 

Note that the current isochrones have been transformed into 
the observational plane usin g the homogeneous set of a tmo- 
sphere models provided by IBrott & Hauschildt I (I2005L see 
Sect. |4]for more details). To constrain th e difference between 
the ma gnit udes and colors base d on the IBrott & Hauschildtl 
(120051) and I Allard et al. I (1201 ll) sets of atmosphere models 
we performed a test. We transformed the best fit cluster 
isochrone (12 Gyr) and chemical compo sition ([Fe/Hl=-0.87 
dex, [c^e]=-(-0.3) using also the recent lAllard et alH (1201 ll) 
spectrjn. 

To overcome problems in the adopted solar mixtures the 
two sets of atmosphere models were interpolated at fixed 
metallicity ([Fe/H]=-0.87 dex). We found that the difference 
in the r' magnitude is negligible, and indeed, it is of the order 
of 1% at MSTO and becomes at most of the order of a few 
hundredths of magnitude in the region across the MSK. We 
performed the same experiment, but the atmosphere models 
were interpolated at fixed global metallicity Z. The difference 
between the two sets of atmosphere models was, once again, 
minimal. Thus suggesting that new atmosphere models have 
a minimal impact on the conclusions of the current investiga¬ 
tion. A more detailed investigation of the the impact that CIA 

(Allard et al. I <20111) spectra are available at the following URL: 
http: //phoenix. ens- lyon. £r/Grids/BT- Settl/ 
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and molecular line opacities have on UV, optical and NIR col¬ 
ors over a broad range of metal abundances will be addressed 
in a forthcoming paper. 

6. THE ABSOLUTE AGE OF M71 

Current photometry is very accurate and precise over at 
least ten magnitudes (see Fig. 5). This outcome applies in 
particular to the bluer bands (u',g',r'). Moreover, theory 
and observations display a well defined knee (r' ~21.5 mag) 
in the r', u' - g' (u' - g' ~2.9 mag) and in the r', g' - / 
(g' - r' ~1.4 mag) CMP. The knee w as also detected in the 
r',g' - r' CMD bv IClem et al.l (l2008l see their Fig. 12) us¬ 
ing the same data set, but a different data reduction strategy. 
The plausibility of the detection was also supported by the 
knee showed by metal-r ich cluster isochro nes (VandenBerg 
et al. 2006) adopted by IClem et ^ (l2008l see their Fig. 6) 
in the r', u' - g' CMD for the typical absolute age (~12 Gyr) 
and metal-abundance ([Fe/H]—0.71 dex) of M71. A similar 
evidence, but based on multiband (SPSS) observations, was 
brought forward for metal-rich clusters by lAn et al. 1 (l2009t 
see their Figures 5 and 8). 

Therefore, we decided to take advantage of this evidence 
and to estimate the absolute age of M71 by using both the 
classical MSTO method and the new method. The lat¬ 
ter method appears very promising, since the vertical method 
might be more prone to possible systematic uncertainties in 
the metal-rich regime. The HB morphology of metal-rich 
GCs is typically characterized by a stub of red stars. These 
stars in several optical and optical-NIR CMDs appear tilted 
(iRaimondo et al. 1120021 iMilone et al. 1120125 iLagioia et al. I 
1201 4 VI3). Moreover, metal-rich GCs-alike M71-harbour 
either no RR Lyrae or at most a few (only one in 47 Tuc). 
If pre sent they ar e quite often eithe r evolved (iBono et al. I 
l2003h or peculiar (iPritzl et al. Il2000ll . The vertical method 
is typically anchored to the mean magnitude of cluster RR 
Lyrae. Therefore, the application to metal-rich clusters is- 
once again-difficult, since a theoretical/empirical correction 
is required. The homogeneity of the diagnostic adopted to es¬ 
timate the absolute age is a crucial issue in this context, since 
we are interested in constraining possible differences in age 
between metal-poor and metal-rich GCs. 

The approach adopted to estimate the observed magnitudes 
of bot h MSTO and MSK has already been discussed in de¬ 
tail bv iBono et al. I (1201 Oh in their investigation of NGC 3201. 
In the following, we briefly outline the main steps of the 
quoted approach. The observed ridgeline of the r' ,g' - r' 
CMD was equally sampled (A^/ =0.01 mag) using a cubic 
spline, and then, the MSTO and the MSK magnitudes were 
estimated as the points showing the minimum curvature along 
the MS ridgeline. We found that the MSTO is located at r' - 
17.79 + 0.01 mag, while the MSK at r' = 20.66 + 0.01 mag. 
Note that the uncertainties account for the photometric errors 
and for errors in the ridgeline and in the location of the above 
points (see Fig.j^i. 

To estimate the absolute age of M71 with the classical 
MSTO method we adopted the more metal-poor (Z=0.003, 
7=0.254) set of isochrones, since they take account for the 
observed stars in the different CMDs. Note that the marginal 
difference in color between observations and theory-in the r', 
u' - g' CMD-has a minimal impact on the absolute age esti¬ 
mate, since we are only using the magnitude of the MSTO. 

Data plotted in our Fig. Vindicate an absolute age for M71 
of order of 12 Gyr. To properly estimate the age and the 
errors affecting the ages based on the MSTO we followed 


the approach suggested bviRenzini & Fusi Peed I (Il988h (see 
also iBuonanno et ainil99^ . We selected the r'-band abso¬ 
lute magnitude of the MSTO on cluster isochrones for ages 
ranging from 8 to 14 Gyr. These estimates were performed 
in three sets of r' ,g' - / isochrones computed assuming dif¬ 
ferent chemical compositions, namely: Z=0.002, 7=0.252; 
Z=0.003, 7=0.254 and Z=0.004, 7=0.256. Note that we 
adopted this CMD, because these are the bands in which both 
theory and observations are more precise. On the basis of 
the quoted estimates we performed a linear regression among 
age, metallicity and r'-band absolute magnitude of MSTO. 
We found the following relation: 

log f[Gyr] = (-1.06 + 0.06) -h (0.51 + 0.01) x M™ - 

(0.09 + 0.01) X [Fe/H] (2) 

where t is the absolute age in Gyr and the other sym¬ 
bols have their usual meaning. Using the above relation 
we estimated an absolute age for M71 of 12+2 Gyr ( see 
top panel of Fig|9] ). The error budget of this determina¬ 
tion takes account of photometric errors in the magnitude 
of the MSTO {o-f ~0.01-0.02 mag), in the cluster red¬ 
dening (cr£(^_yj ~0.05 mag), in the true distance mod¬ 
ulus {cTfi ~ 0.15 mag) and in the chemical composition 
(o‘[M//r] ~0.1 dex). The current estimate is, within the er¬ 
rors, in good agreement with the absolute ages for M71 
available in the literature (see Table |2|. The agreement is 
quite good with recent estimates based both on ground-based 
dGrundahl, Stetson & Andersen 11200^ iBrasseur et al. l[20T^ 

and on space (VI 3) photometry together with up-to-dated 
cluster isochrones (I VandenBerg et al.l2012i) . 

To estimate the absolute age with the method of the 
we adopted the same sets of cluster isochrones (Z=0.003) we 
used for the MSTO method. However, the cluster isochrones 
with ages ranging from 8 to 14 Gyr were equally sampled 
using a cubic spline. The same algorithm adopted to esti¬ 
mate the observed MSTO and the MSK on the ridgelines was 
also adopted on the splined isochrones. The observed and 
predicted differences in magnitude (A^™^^) between the MSK 
and the MSTO (r', r' - g') are plotted in the bottom panel of 
Fig.|9l The predicted A;^“®, at fixed metal content, plotted 
in this figure show a well defined linear trend over the inves¬ 
tigated age range. To provide solid constraints on the errors 
affecting the absolute age based on the we followed the 

same approach adopted for the MSTO. Using the r',g' - r' 
cluster isochrones, we obtained the following relation: 

log f[Gyr] = (2.18 + 0.06) - (0.48 + 0.02) x A^^" - 

(0.30 ± 0.02) X [Fe/H] (3) 

where the symbols have their usual meaning. Using the above 
analytical relation, we found an absolute age for M71 of 
12 Gyr. The error on the current age is +1 Gyr and takes ac¬ 
count of uncertainties in photometry and in iron abundance. 
The above error is roughly a factor of two smaller than the 
typical error of absolute ages based on the MSTO. The differ¬ 
ence is due to the fact that the A^^^® diagnostic is independent 
of uncertainties affecting both the distance modulus and the 
cluster reddening. Moreover, the difference in magnitude- 
A^“®-is also minimally affected by uncertainties in the abso¬ 
lute photometric zero-points. This means that this approach 
can also provide very accurate estimates of the relative ages 
of stellar systems. 
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In this context, it is worth mentioning that the coefficient of 
the age indicator attains, as expected, very similar values be¬ 
tween MSTO and This is the consequence of the fact 

that the above diagnostics rely on the same age indicator: the 
MSTO. On the other hand, the dependence on the iron abun¬ 
dance is larger for the than for the MSTO (0.3 vs 0.1 
dex). The difference is independent of the adopted magnitude 
(r') and color {g' - r'). We performed several tests using the 
same cluster isochrones and we found that the metallicity de¬ 
pendence is strongly correlated to the definition of the MSK 
(maximum vs minimum curvature). Moreover, we also found 
that the metallicity dependence of the age diagnostic (A!^q“) 
becomes similar to the MSTO using the the difference in color 
instead of the difference in magnitude. The above evidence 
suggests that the precision of the different age diagnostics is 
a complex balance among photometric and spectroscopic pre¬ 
cision age sensitivity together with uncertainties affecting dis¬ 
tances moduli and reddening corrections. We plan to provide 
a detailed mapping of the sensitivity of optical, optical-NIR 
and NIR CMDs in a forthcoming investigation. 

7. SUMMARY AND FINAL REMARKS 

We performed accurate and deep multiband r', i',z') 

photometry of the metal-rich Galactic globular M71 
(NGC 6838). The images were collected with MegaCam at 
CFHT and cover one square degree around the center of the 
cluster. The image quality and the strategy adopted to perform 
the photometry on both shallow and deep exposures allowed 
us to provide precise photometry in all the bands from the tip 
of the RGB to five magnitudes fainter than the MSTO. The 
extend ed data set of local standards provided by IClem et aH 
(l2007h allowed us to fix the absolute zero-point of the above 
bands with a precision that is, on average, better than 0.02 
mag. 

The selected cluster is projected onto the Galactic bulge. 
This means that absolute age estimates are hampered by 
the contamination of field bulge stars. To overcome this 
thorny problem we deviced a new approach based on the 
(r',u' - g',g' - r') Color-Color-Magnitude Diagram. First we 
derived accurate ridgelines of candidate cluster stars in the 
r',u'-g' and in the r',g'-r' CMDs using iso-density contours. 
Then, they were combined to provide a 3D trace to properly 
selected candidate cluster stars. The key advantage of this ap¬ 
proach is that it takes advantage of both intrinsic properties 
(color selections) and of the clustering (same distance). The 
results appear very promising, but the approach can be further 
improved using optical-NIR color selections. 

We computed specific sets of cluster isochrones covering a 
broad range in ages and in the adopted chemical composition. 
Moreover, to properly fit observed red giant branch stars, the 
evolutionary models were computed by adopting three dif¬ 
ferent values for the mixing length parameter. To constrain 
possible systematics in the transformation into the observa¬ 
tional plane, we paid special attention on the adopted stellar 
atmosphere models and in the adopted-pass bands. We found 
that the latter plays a crucial role in the comparison between 
theory and observations. 

We also investigated the impact that the collisional induced 
absorption (CIA) opacity has in occurrence of the main se¬ 
quence knee in the Sloan bands. We found that CIA affects not 
only NIR magnitudes, but also blue bands such as the u' and 
the g' band. The sensitivity of the short wavelength regim e on 
the CIA opacity was predicted bv iBorvsow et al. I (119971) and 
now soundly confirmed by both theory and observations. 


The comparison between theory and observation is quite 
good over the entire magnitude range covered by observations 
in the different CMDs. There is evidence of a mild difference- 
a few hundredths of magnitude-between cluster isochrones 
and observations. However, it is not clear whether the differ¬ 
ence is due either to the adopted atmosphere models or to the 
adopted reddening law. 

We estimated the absolute age of M71 using two different 
diagnostics. Using the MSTO and the cluster isochrones in 
the r',g' - r' CMD, we found an age of 12+2 Gyr. The er¬ 
ror takes account of uncertainties in cluster distance, redden¬ 
ing and iron abundances together with photometric uncertain- 
ties. We performed the same estimate using the knee-method 
(iBono et al. 11201 Oh . i.e. the difference in magnitude between 
the MSTO and the MSK, in the same CMD. We found a clus¬ 
ter age that is identical to the age based on the MSTO, but 
the uncertainty is of a factor of two smaller. The difference is 
mainly due to the fact that the adopted diagnostic—A^ q“— is 
independent of uncertainties in cluster distance and redden¬ 
ing. However, it is more prone to uncertainties affecting clus¬ 
ter iron abundance. 

The above age estimates support, within the errors, re¬ 
cent age estimates of metal-rich Galactic globulars. One 
of the most m etal-rich Galac t ic glo bular-NGC 6528-was 
investigated by iLagioia et al. I (1201 4l) using deep ACS and 
WFC3 images. Using the MSTO and a similar set of clus¬ 
ter isochrones they found and absolute age of 11 + 1 Gyr. 
The absol ute age of this cluster was independently con- 
hrmed by ICalamida et al. I (1201 4ft using for the same clus¬ 
ter accurate ground-based Stroemegren and NIR photome¬ 
try (t=ll + l Gyr) . The outcome applies to the age esti¬ 
mates of met al-ri ch bulge clusters provided bvIZoccali et al. I 
(l2003l l2004ft. bv[Dotter. Saraiedini & Anderson! (1201 Ih and 
bv lBellinVeTal.feoTl . 

The scenario emerging from the above empirical evidence 
is that metal-rich Galactic globulars appear to be, within 
the errors, coeval with meta l-poor globulars (iMonelli et al. I 
1201 3t iDi Cecco et al. I l2010ft . This means that the quoted 
globulars do not show a clear evidence of an age metallic¬ 
ity relation. This is a preliminary conclusion, since we still 
lack homogeneous age estimates of a sizeable sample metal- 
rich clusters with a precision better than 1 Gyr. The fact that a 
significant fraction of metal-rich clusters are located inside or 
projected onto the Bulge is calling for new precise and deep 
NIR photometry. Together with the obvious bonus concerning 
the reduced impact of uncertainties in the cluster reddenings 
and the possible occurrence of differential reddening, there is 
also the advantage to fully exploit the precision of the new 
diagnostic (A;^*^®). 

Finally, we would like to stress that the above findings ap¬ 
pear also very promising concerning the solid evidence of 
a dichotomic age distribution of Galactic globulars . This 
evidence was brought forward by ISalaris & Weiss I (l2002h 
and su bsequently soundly co nhrmed bv iMar in-Franc h et al. I 
(l2009ft an more recently, by [VandenBerg et al. I (120131) and 
iLeaman, VandenBerg & Mendel 1(120131) . The use of the 
parameter appears even more promising in this context ap¬ 
pears even more promising. The reason is twofold, a) em¬ 
pirical evidence indicate that the morphology of the MSK is 
well defined when moving from the meta l-rich (NGC 6528 
ISaraiedini et al. Il2007t ILagioia et al. 1120 in preparation) to 
the metal-intermedi ate (Bono et al. 2010 ) and to the metal- 
poor regime (Ml5. iMonelli et £171120151) . b) Relative ages 
based on the are less prone to systematics, since the 
















































8 


Di Cecco et al. 


MSK is independent of cluster age and mildly affected by 
changes in helium content. The vertical method that is based 
on the difference between the HB luminosity level and the 
MSTO is more prone to the quoted possible systematics. In¬ 
deed, the HB luminosity level is affected by both cluster age 
and helium content. 

The suggested experiment appears to be luckily supported 
by the unique opportunity to use in the near future adap¬ 
tive optics systems at the 8 m class telescope with superb 
image quality and spa tial resolution (ISchreiber et al. Il2014t 
iFiorentino et al. I^2014^ . 
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TABLE 1 

Ridgelines of M71 in the r', u' — g’ and r', g' — r' CMDs“. 


r' 

mag 

0-r> 

mag 

u' ~g' 
mag 

mag 

g' -r' 
mag 

mag 

11.49 

0.01 

3.94 

0.04 

1.58 

0.12 

11.69 

0.01 

3.70 

0.04 

1.47 

0.04 

11.89 

0.01 

3.48 

0.07 

1.38 

0.09 

12.09 

0.01 

3.30 

0.06 

1.32 

0.10 

12.21 

0.01 

3.21 

0.06 

1.28 

0.06 

12.33 

0.01 

3.13 

0.04 

1.25 

0.04 

12.50 

0.01 

3.02 

0.03 

1.21 

0.03 

12.66 

0.01 

2.92 

0.08 

1.18 

0.05 

12.82 

0.01 

2.83 

0.06 

1.15 

0.04 

13.02 

0.01 

2.72 

0.03 

1.11 

0.03 

13.22 

0.01 

2.62 

0.04 

1.08 

0.03 

13.46 

0.01 

2.52 

0.03 

1.05 

0.04 

13.70 

0.01 

2.42 

0.03 

1.01 

0.03 

13.94 

0.01 

2.32 

0.03 

0.98 

0.02 

14.18 

0.01 

2.24 

0.02 
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0.02 

14.42 

0.01 
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0.02 

14.67 

0.01 
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0.02 
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14.91 

0.01 
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0.01 

1.98 

0.02 
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0.02 
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0.01 
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0.02 

0.86 

0.02 

15.63 

0.01 
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0.02 
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0.02 

15.91 

0.01 
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0.02 

0.83 

0.02 
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0.01 
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0.02 

0.82 

0.02 

16.19 

0.01 
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0.02 

0.82 

0.02 
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0.01 

1.80 

0.01 
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0.01 
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0.01 
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0.01 
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1.76 
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0.01 
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0.62 
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0.61 

0.01 
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1.25 
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TABLE 1 — Continued 


r' 

mag 

(Tr’ 

mag 

mag 

mag 

mag 

(Tgt-f 

mag 

17.57 

0.01 

1.24 

0.01 

0.60 

0.01 

17.59 

0.01 

1.24 

0.01 

0.60 

0.01 

17.61 

0.01 

1.24 

0.01 

0.60 

0.01 

17.63 

0.01 

1.24 

0.01 

0.60 

0.01 

17.65 

0.01 

1.23 

0.01 

0.60 

0.01 

17.67 

0.01 

1.23 

0.01 

0.60 

0.01 

17.71 

0.01 

1.23 

0.01 

0.60 

0.01 

17.75 

0.01 

1.23 

0.01 

0.60 

0.01 

17.79 

0.01 

1.23 

0.01 

0.60 

0.01 

17.83 

0.01 

1.23 

0.01 

0.60 

0.01 

17.89 

0.01 

1.23 

0.01 

0.60 

0.01 

17.95 

0.01 

1.24 

0.01 

0.60 

0.01 

18.07 

0.01 

1.25 

0.01 

0.60 

0.01 

18.23 

0.01 

1.27 

0.01 

0.61 

0.01 

18.39 

0.01 

1.30 

0.01 

0.62 

0.01 

18.55 

0.01 

1.33 

0.01 

0.64 

0.01 

18.71 

0.01 

1.38 

0.01 

0.66 

0.01 

18.87 

0.01 

1.43 

0.01 

0.68 

0.01 

19.03 

0.01 

1.49 

0.01 

0.71 

0.01 

19.23 

0.01 

1.57 

0.01 

0.74 

0.01 

19.44 

0.01 

1.67 

0.01 

0.78 

0.01 

19.64 

0.01 

1.78 

0.01 

0.83 

0.01 

19.84 

0.01 

1.90 

0.01 

0.88 

0.01 

20.04 

0.01 

2.03 

0.01 

0.93 

0.01 

20.20 

0.01 

2.14 

0.01 

0.97 

0.01 

20.40 

0.01 

2.29 

0.01 

1.04 

0.01 

20.60 

0.01 

2.43 

0.01 

1.11 

0.01 

20.88 

0.01 

2.59 

0.01 

1.20 

0.01 

21.12 

0.01 

2.73 

0.01 

1.28 

0.01 

21.36 

0.01 

2.87 

0.01 

1.35 

0.01 

21.56 

0.02 

2.97 

0.01 

1.41 

0.01 

21.76 

0.02 

3.07 

0.02 

1.46 

0.02 

21.84 

0.02 

3.11 

0.02 

1.48 

0.02 


^ The sampling in magnitude and in color along the r' ,u' — g' and the r',g' — r' ridgelines is not uniform. The sampling increases in the regions in which 
there are relevant changes in the slope (main sequence turn off, base of the red giant). The cr marks the uncertainty either in magnitude or in color (summed in 
quadrature). 


TABLE 2 

Absolute ages, distance moduli, reddening estimates and iron abundances for M71. 


DMv“ 

E{B-Vf 

Age‘= 

[Fe/H]"* 

Notes'^ 

mag 

mag 

Gyr 



13.70 

0.28 

14±2; 16±2 

-0.78 

H92 

14.09±0.15 

0.28 

8±1 

-0.70 

R92 

13.60±0.10 

0.27±0.05 

18 

-1.02 

GOO 

13.71±0.04±0.1 

0.28 

12 

-0.70 

G02 

13.78 

0.20 

11 

-0.80 

BIO 

13.69 

0.24 

11.00 ±0.38 

-0.82 

V13 

13.84 

0.25 

12±2 

-0.78 

MSTO 



12±1 

-0.78 

A Knee 
^TO 


^ Apparent distance modulus and its error when estimated by the quoted authors. 

^ Adopted cluster reddening. 

^ Absolute cluster age. 

^ Adopted iron abundance. 

^ Notes: H92, IHodder'etaTI fT9^ : R92. [Mdl fT99^ : GOO. IGeffert & Maintz I j2QOOi) : G02, IGrundahl. Stetson & Andersen I <200l ; BlQ. IBrasseur et all 
nnwi) : V13. IVandenBerg et al. M2Q13I) . MSTO: age estimate based on the main sequence turn off; age estimate based on the difference in magnitude 

between MSTO and knee. 


TABLE 3 

Effective temperatures and gravities of the selected models to study the effect of CIA. 


Tcfs 

logg 

[K] 

[cgs] 

3777 

4.98 

3904 

4.91 

4081 

4.82 
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TABLE 4 — Continued 
CIA Ref. 


TABLE 4 

CoLLisioNAL Induced Absorption opacities included in the PHOENIX BT-Settl 
ATMOSPHERE MODELS (aLLARD ET AL. 2011). 


CIA 

Ref. 

H 2 -H 2 

1,2 

H 2 -He 

3,4 

H 2 -H 

5 

He-H 

6 

H 2 -CH 4 

7, 8 

H 2 -N 2 

9, 10 

N 2 -CH 4 

11 , 12 

N 2 -N 2 

13,14 

CH 4 -CH 4 

15 

CO 2 -CO 2 

16, 17 

H 2 -Ar 

18, 19 

CH 4 -Ar 

20,21 


* IBorvsowl 420021) : ^ lAbel et al~l 420111) : ^ IBorvsow. Frommhold & Moraldi I 419891) : ^ IBorvsow et al. I 419971) : ^ IGustafsson & Erommhold I 420031) : 
® IGustafsson & Frommhold I 420011) : ^ IBorvsow & Frommhold I 4198614) : ® IBorvsow. Frommhold & Dorel 419861) : ® IBorvsow & Frommhold I 419861) : 
IDore. Borvsow & FrommhoCTI 419861) : ** IBorvsow & Tan^ 419931) : IBimhaum. Borvsow. & Buechele I 419931) : IBorvsow & Frommholdl 41986cl) : 
IBorvsow & Frommhold I fTW7^ : ‘^ IBorvsow & Frommholdl41987bl) : *^ IGruszka & Borvsow I41997^ : '^ IGniszka & Borvsow 141993) : *^ IMever & Frommholdl 
419861) : '^ IBorvsow & Moraldi 1 419941) : IDore & Filahozzn4199(ll) : ^* IBorvsow & Moraldi 1 4 1 993l) . 
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u'-g' g'“'' 9'“Z' a^. 


Fig. 1.— From left to right: four CMDs of M7 1 based on optical images collected with the MegaCam at CFHT. Data plotted in the different CMDs were selected 
using photometri c quality paramete rs « 1, \sharp\ < 1, (Tg'-/ < 0.07) and according to the radial distance 60" < r < 250" . The red lines display the fiducial 
lines provided bv IClem et alj 420081) . while the red bars plotted on the left side take account of the photometric error both in magnitude and in color. The sample 
of stars located at r'~14.2, u' - g'~\.15 mag are cluster red HB stars. Data plotted in the rightmost panel display the photometric error in the r'-band. 
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1-5 2.5 25 0,5 1-0 1.5 20 


u'-g' g'-r‘ 


Fig. 2.— Left - Iso-density contours (black lines) of the entire data set in the r',u' — g' CMD. The ridgeline of the candidate cluster stars is overplotted with red 
asterisks. Right — same as the left, but for the r',g' - r' CMD. 
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Fig. 3.— r',u' — g' ,g' — r' Color-Color-Magnitude-Diagram (CCMD) of M71. The blue dots are the candidate cluster stai's selected along the 3D ridgeline (red 
line). The red dots are candidate field stars. 
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1 2 


1 4 - 


1 6 


1 8 - 


20 


22 


Totq I 

1 0''<r<450’ 

CTg<0.05 

a^<0.05 . 



M71 

S e I e c t e d 
H B 



Field 





0.5 1.0 1.5 

g'-r' 


0.5 1.0 1.5 

g'-r’ 


Fig. 4.— Left — r',g' - r' CMD of the entire photometric catalog. Stars plotted in this CMD were selected in photometric error (rr^ < 0.05, cr^ < 0.05) and in 
cluster radial distance (10" < r < 450" ). Middle - Same as the left, but for candidate cluster stars. The blue dots display cluster red HB stars. Right - Same as 
the left, but for candidate field stars. 
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Fig. 5.— From left to right comparison between theory and observations in different CMDs of M71. Black dots display candidate cluster stars selected in 
photometric en'or and in cluster radial distance. The red lines display two cr-enhanced cluster isochrones at fixed metal content (Z=0.003, Y=0.254) and ages of 
11 and 13 Gyr, respectively. The tme distance modulus and the cluster reddening adopted to overplot evolutionary prescriptions are also labelled. The reddening 
in the different bands was estimated using the semi-empirical Cardelli’s relation. The black arrows show the reddening vectors in the different CMDs for an 
arbitrary change in cluster reddening. 
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ij'—g' g’ — r' g’ — i’ g’ — z' 

Fig. 6.— Same as Fig. 5, but the comparison is with two more metal-rich (Z=0.004, Y=0.256) a-enhanced cluster isochrones. 



Fig. 7.— Comparison between cluster isochrone (12 Gyr) transformed into the observational plane adopting stellar atmosphere models either including (red 
line) or neglecting (black line) CIA opacity. 
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Fig. 8. — From left to right optical CMDs of candidate cluster stars. Stars plotted in the dilferent CMDs were selected according to photometric error and cluster 
radial distance. The red lines display the ridgelines, while the two large blue filled circles plotted in the r',g' — r' CMD mark the position of both the MSTO and 
of the MSK. 
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Fig. 9. — Top - Predicted absolute magnitude of the MSTO as a function of cluster age. The MSTO was selected in / ,g' — r' isochrones constructed by adopting 
three different cr-enhanced chemical compositions (see labeled values). The observed unreddened MSTO value is 4.07 mag and it is shown (by the vertical 
bar) considering the uncertainties of the photometiy and of the distance modulus and reddening. The horizontal bar shows the estimated uncertainty on the age. 
Bottom - Same as the top, but for the pai'ameter (2.87 mag). The vertical bar only takes account of the photometi'ic error. The horizontal bar shows the 

estimated uncertainty on the age. 








